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ABSTRACT 

We study the coUisional evolution of km-sized planetesimals in tight binary star sys- 
tems to investigate whether accretion towards protoplanets can proceed despite the 
strong gravitational perturbations from the secondary star. The orbits of planetesimals 
are numerically integrated in two dimensions under the influence of the two stars and 
gas drag. The masses and orbits of the planetesimals are allowed to evolve due to colli- 
sions with other planetesimals and accretion of collisional debris. In addition, the mass 
in debris can evolve due to planetesimal-planetesimal collisions and the creation of new 
planetesimals. We show that it is possible in principle for km-sized planetesimals to 
grow by two orders of magnitude in size if the efficiency of planetesimal formation is 
relatively low. We discuss the limitations of our two-dimensional approach. 
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1 INTRODUCTION 

We consider the problem of planetesimal evolution in tight 
binary star systems such as 7 Cephei and a Centauri. To 
date several hundred extrasolar planets have been detected 
- 20% of which orbit the primary of a binary or multiple 
system ( [Desidera fc BarbierilpOOT] ). In most of these cases 
the stars are widely separated and do not significantly af- 
fect the evolution of planets, but 7 Cephei, G186, HD41004A, 
and possibly HD196885 are examples of tight binary systems 
with a semi-major axis, Ob, around 20 AU, that have Jupiter 
mass planets. 7 Cephei has the most extreme binary param- 
eters of the known tight binaries with detected planets. The 
primary star is a K giant (Ma* ~ 1.4 M©) the secondary is 
a M dwarf (A/b* ~ 0.4 M0), the binary has a high eccen- 
tricity (eb = 0.4), and a small Ob of 20 AU. The planet is 
1.6 Mj with an a ~ 2 AU ( [Neuhauser et al.||2007| . Due to 
its extreme characteristics 7 Cephei is often considered one 
of the most stringent tests for planet formation models. 



Previous work by Thebault et al. ( 2004 1 and Quintana 



et al. (20071 has shown that planets can form in 7 Cephei 



if planetary embryos can form within the stable region < 3 
AU ( Wiegert fc Holman|1997 l. However, work on the earlier 
stage of planet formation has shown that it is very difficult 
to grow planetesimals that are initially in the km size regime 
because of the secular perturbations from the secondary star 
on the planetesimals and the gas disc. The coupled effect of 
perturbations from the binary and the gas drag from the 
disc causes differential orbital phasing, and high speed de- 
structive impacts between planetesimals of different sizes. 
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Thebault et al. ( 2006 1 suggest that the destructive collisions 
prohibit the growth of planetesimals in 7 Cephei at 2 AU, 
the current location of the planet. 

Although planets have yet to be detected in the 
a Centauri system, the proximity of the system and the 
roughly solar mass of both stellar components make the 
system a prime target for finding an Earth mass planet 
in the habitable zone (Guedes et al. 20081. However, the 
a Centauri system is more extreme than 7 Cephei with a 
mass ratio close to one (Ma* = 1.1 M© and Ms* = 0.93 
Mq), eb = 0.52, and Ob = 23.4 AU. Thebault et al.| (|2008 



|2009[ ) found that impact speeds between different sized plan- 
etesimals would be above the disruption threshold due to dif- 
ferential orbital phasing, and planetesimals could not grow 
in the habitable zone in either a Centauri A or B. 

In this letter we show that planetesimal growth from 1 
to 100 km is possible in principle in these tight binaries when 
the collisional evolution of the planetesimal disc is taken into 
account. Disruptive collisions will prevent differential orbital 
phasing from being set up and, at the same time, create a 
reservoir of small debris that can be accreted onto the left 
over planetesimals. 



2 DIFFERENTIAL ORBITAL PHASING 

In the absence of any gas, a tight eccentric binary companion 
will excite eccentricities of particles orbiting the primary. 
For particles starting on circular orbits, their eccentricity 
will oscillate around the forced eccentricity a: 
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Figure 1. Eccentricity vs. semi-major axis of coUisionless, mass- 
less test particles in the 7 Cephei system after 100 binary orbits 
in a gas disc with a constant gas density of pg = 1.4 • 10~^ g 
cm~'^. Colour indicates size in km; the black curve with the large 
oscillations denotes a gas-free case. 



sion. The latter will be of the order of the escape velocity, 
VcBc = ^JlGM I R. Focusing on a = 1 AU around a So- 
lar type star we have Wesc = 4.5 • 10~^(i?/km)an, where Q, 
is the angular velocity. The number density n is given by 
n = Es/(AzM), where Es is the surface density of solids 
and Az is the vertical extent of the particle disc. We can 
write Az ~ ai, with i the maximum inclination, which we 
can assume to be small. In fact, we expect 2j « e ~ Vas^jaQ.. 
Plugging this all in we find for the collision time scale: 
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where pp is the bulk density of the particles. This means that 
for typical parameters at 1 AU, a km-sized planetesimal will 
undergo a collision once every few thousand years. 



with amplitude 2ef, where a is the semi-major axis of the 
particle ( Heppenheimer|[l978 1 . Since there is no dissipative 
force in the system, the oscillations do not damp and their 
spatial frequency increases with time, which eventually leads 



to orbital crossing of neighbouring particles ( Thebault et al. 
2006 1 . This is illustrated by the black curve in Fig.[lJ As soon 
as particles at an eccentricity maximum can collide with 
particles at an eccentricity minimum the resulting encounter 
velocities are too high for accretion to occur. 

The presence of a gas disc and the associated aerody- 
namic gas drag on the particles can effectively damp the os- 
cillations leading to a well-defined equilibrium eccentricity 
distribution eo — eo(a) (see [Paardekooper et al.|2008[ ), with 
corresponding periastron alignment. The time scale for this 
equilibrium to be reached is the secular time scale. There- 
fore, it takes many binary orbits before planetesimals are 
on their equilibrium orbits. Once they have reached their 
equilibrium orbits particles of the same size have orbits that 
are in phase and, therefore, they will not suffer from high- 



speed collisions with each other (Marzari & SchoU 20001 



However, since the magnitude of the gas drag force depends 
on the particle size, each particle size will have a different 
eo(a), a phenomenon that is called differential orbital phas- 
ing ( Thebault et al.|2006 l. In Fig. [l] particles of a single size 
(and colour) follow phased orbits with no oscillations. The 
residual oscillations that can be seen for the largest particles 
will eventually be damped completely by gas drag. However, 
for particles of differing size there is a considerable spread 
in eccentricity and, not shown, longitude of periastron. The 
difference in eccentricity magnitude can easily exceed 0.01 
for particles of 1 and 5 km. This implies that particles of dif- 
ferent size will undergo high-speed collisions, which are usu- 
ally destructive. As a result, it is difficult to form planetary 
cores in tight binary systems if the only growth mechanism 
is accretionary planetesimal-planetesimal collisions. 



3 COLLISION TIME SCALE 
3.1 Single star case 

Consider a population of equal-sized bodies of mass M 
around a single star at 1 AU. The collision time scale is 
given by Tc — 1/ (nR^n/^v), where R is the radius of the 
bodies, n the number density and An the velocity disper- 



3.2 Binary star case 

Perturbations due to a coplanar, eccentric binary will in- 
crease the eccentricity of the planetesimals. Differential or- 
bital phasing can easily give rise to the velocity dispersion 
entering the collision time scale will go up by two orders 
of magnitude compared to the single star case. Keeping all 
other parameters the same, a km-sized planetesimal will now 
undergo a collision every 10-20 orbits at 1 AU. If we take the 
binary companion to have a semi-major axis of 20 AU, this 
collision time scale amounts to 25 % of a binary orbit. This 
should be compared to the time scale for a particle to reach 
its equilibrium orbit in the binary system, which happens on 
a secular time scale (many binary orbits). Therefore, phys- 
ical collisions are of crucial importance for the evolution of 
the planetesimal population. It is important to stress that 
not only, as has been realised before, are collisions destruc- 
tive when differential orbital phasing happens, they can ac- 
tually prevent this size-dependent orbital structure from be- 
ing set-up in the first place. It is therefore necessary to take 
collisions into account when studying the effect of differen- 
tial orbital phasing on planetesimal accretion. 

Another consequence of having many possibly catas- 
trophic collisions in the system is a large amount of small 
debris. This small debris can act as a source of new planetes- 
imals or be accreted onto remaining larger bodies. All these 
ingredients need to be taken into account in a consistent 
model of planetesimal accretion in binary systems. 



3.3 Two-dimensional disc 

In order to reduce the computational cost of the models, we 
will work in a two-dimensional (2D) geometry. The collision 
time scale, for orbits with random phase, is then given by 
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where n is the surface number density and e is a measure 
of the velocity dispersion (e ^ 4.5 • lO"'^ ff Av = Vesc for 1 
km sized bodies). The collision time scale for a 2D disc is 
artificially short; we have to correct for this by choosing a 
low surface density in the simulations so that for the unper- 
turbed disc, Tc 2D ~ Tc- 
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4 MODEL DESIGN 



4.4 Collision outcomes 



From the above discussion it is clear that the effect of coUi- 
sions, on the size distribution as well as on the orbital ele- 
ments of the planetesimals, can not be ignored. As a result, 
we use a simple model that incorporates the three necessary 
components: gas, planetesimals and small dust. 



4.1 Gas disc 

For simplicity, we take the gas disc to be static and circular. 
Although the gas disc is expected to become eccentric under 



We use the velocity-dependent catastrophic disruption cri- 



the influence of the binary companion ( Paardekooper et al 



2008 Kley fc Nelson 20081 the qualitative outcome of the 



model presented here is not affected by the assumption of 
zero eccentricity. 'Paardekoo per et al.| ( |2008[ ) showed that if 
the gas disc does not follow the forced eccentricity, differ- 
ential orbital phasing will occur, which always occured in 
the full hydrodynamical simulations. A circular gas disc is, 
therefore, a reasonable starting point. 

The gas density used in this work is assumed to be con- 
stant with a. The value is chosen to be close to the Minimum 



Mass Solar Nebula at 1 AU, 
resulting gas drag force is 



pi = 1.4-10 g cm 



-fdrag — 
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with Cd the drag coefficient. We take Cd = 0.4, appropriate 
for spherical bodies. 



4.2 Small dust 

We embed the gas disc with a population of tightly coupled 
small particles, which are assumed to move on circular or- 
bits together with the gas. They can be accreted onto exist- 
ing planetesimals, but also form new planetesimals, and are 
created in destructive collisions between planetesimals. The 
small dust particles are distributed over typically 32 radial 
bins. At the start of the simulation, the mass in dust is equal 
to Md = MminN, where Mmin is the mass of a planetesimal 
of the minimum size we consider (usually 1 km), and N is an 
input parameter specifying the total mass. Note that A*' is 
also the maximum number of particles we could have at any 
one time in the simulation. The dust mass inside a radial 
bin is smeared out to give a smooth surface density Ed. 



4.3 Planetesimals 

We model planetesimals as test particles, moving under the 
influence of gravity from both stars as well as gas drag. 
Planetesimals have a minimum size of 1 km (everything 
smaller is taken to be small dust), and can grow by low- 
velocity collisions with each other and by accreting small 
dust. Collisions are detected by assigning an inflated radius 
-Rinfl cx i? to the particles and testing whether two particles 
overlap. The planetesimal disc is characterised by the prod- 
uct RinfiN, which, together with Av, sets the collision time 
scale. Since Av is set by the gravitational perturbations due 
to the binary companion and is independent of A'^ and Tiinfl, 
we expect simulations with the same value for Ri^RN to give 
similar results, a well-known result of collision models. 



terium of Stewart & Leinhardt ( 2009 1 (strong particle ver- 
sion) to determine the size of the largest remnant when a 
collision is detected. We then use the technique described in 



Wyatt & Dent ( 2002 1 to predict the second-largest remnant. 



The total mass minus the largest remnant is assumed to fol- 
low a power law size distribution with index —1.93, from 
which the number of bodies larger than R, N{> R), can be 
derived. The size for which N{> R) — 2 is the size we take 
for the second largest remnant ( Wyatt fc Dent||2002 l. Any 
mass that is below the minimum particle size is added to the 
small dust. To keep the total number of particles tractable, 
we do not keep the third largest remnant, even if it is bigger 
than the minimum size. 



4.5 Planetesimal formation 

Contrary to all previous studies, we do not let all planetesi- 
mals appear at t = 0, but instead have them form from the 
small dust present in the system. The efficiency of planetes- 
imal formation is a major unknown in this model. We take 
a very simple approach and say that: 



dMA 
dt 
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where tp is an efficiency factor, i?infl,min is the inflated ra- 
dius of the minimum planetesimal size, and Ma is the mass 
available (in a certain radial dust bin) to form planetesimals. 
As soon as Ma > Mmin, we create a new planetesimal. Note 
that for fixed RmnN, Ma oc N, so that the number of plan- 
etesimals at any time is proportional to A''. Increasing A'' 
while keeping RindN constant effectively increases the 'res- 
olution' of the simulation. 

For equal masses in km-sized planetesimals and small 
dust, the collision time scale and the planetesimal formation 
time scale are related through Tp = r^eje-p. If we consider 
the single star case, with e « 10~^ and Tc ~ yr, and if 
we expect planetesimal formation to proceed on a time scale 
of Tp > 10"* yr, we must have tp ~ 10~^. This is the value 
used for the standard model discussed below. 



4.6 Dust accretion 

Planetesimals can accrete dust at a rate 
dM 



dt 



= 2ed-RinflEdju - 



(6) 



where Ed is the surface density of small dust. An efficiency 
factor ed accounts for the possibility that not all dust created 
in collisions is available for accretion, and that the efficiency 
of accreting small bodies may be smaller than 1. Note that 
Ed oc N, so that the amount of dust accreted per planetes- 
imal is the same if N is changed but the product -RinflA is 
kept constant. 



5 RESULTS 

We start by describing a simulation that has N — 10^, 



Rind/ab = 10"^7?/km 
= 1.4 ■ 10 



10" 



Ed 



1 and constant 



® g cm Note that gas drag, and therefore 



4 S.-J. Paardekooper and Z.M. Leinhardt 




Figure 2. Eccentricity vs. semi-major axis of test particles in the 
7 Cephei system after 1000 binary orbits. Color indicates size in 
km. 
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Figure 3. Evolution of maximum particle size (solid curve, top 
panel), mean particle size (dashed curve, top panel), planetesimal 
mass fraction (solid curve, bottom panel) and dust mass fraction 
(dashed curve, bottom panel) for the 7 Cephei system. 



differential orbital phasing, is very strong throughout the 
disc. The minimum particle size is 1 km and there are no 
planetesimals at t — 0. Binary parameters are those of the 
7 Cephei system. Figure [2] shows the resulting (a, e) distri- 
bution after 1000 binary orbits (~ 90000 yr). It is clear that 
significant accretion has taken place in the inner parts of the 
disc, where particles have grown from 1 km to 50 km. Out- 
side a/ab = 0.1, which corresponds to 2 AU, perturbations 
due to the binary are too strong for accretion to occur. 

Inside 2 AU, planetesimals grow to sizes up to 70 km. 
The top panel of Fig.|3]shows the time evolution of the max- 
imum and mean planetesimal size over the whole disc. The 
mean size is dominated by the large number of small plan- 
etesimals in the outer disc, and stays between 10 and 20 km. 
The maximum particle size goes up rapidly by dust accre- 
tion (smooth parts of the curve) and by accreting collisions 
(jumps). The largest planetesimal is destroyed a few times 
as well, but in general the trend is to grow to larger sizes. 
After ~ 600 binary orbits, there is no source of small dust 
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Figure 4. Eccentricity vs. semi-major axis of test particles in the 
a Centauri system after 1000 binary orbits. Color indicates size 
in km. Top panel: a Centauri A. Bottom panel: a Centauri B. 



remaining from which to create new planetesimals, and the 
collision time scale goes up. 

The maximum size that can be reached in this scenario 
is limited by the total amount of solid material present in the 
disc. In the 2D approximation, this mass has to be artificially 
low in order to end up with a realistic collision time scale. 
For a disc that has twice the solid material, keeping all other 
parameters the same, growth up to 100 km was observed. 
Increasing A'^ by a factor of 2 while decreasing i?infl by the 
same factor, which amounts to increasing the resolution of 
the simulation, showed growth up to 150 km. In these higher 
resolution runs more accreting collisions are observed than 
depicted in Fig. |3] This is because the system goes through 
phases of low particle number density, for which the collision 
statistics in the lower resolution runs are not optimal. In this 
sense, Fig. [3] represents a worst-case scenario, where many 
low-velocity collisions are missed, and adding particles, while 
keeping the collision time scale the same, will only favour 
planetesimal growth more. 

Gas drag appears to play only a minor role in deter- 
mining the qualitative outcome of the model. This is mainly 
due to the fact that the planetesimals are weak enough so 
that any small eccentricity difference of ~ 0.01 will lead to 
destructive collisions. Whether this difference is due to dif- 
ferential orbital phasing (under the influence of gas drag) or 
simply orbital crossing (in the absence of gas) does not mat- 
ter. A simulation without any gas showed the same trend as 
in Fig. [3] with growth up to 80 km. 

Crucial parameters are the efficiency of planetesimal 
formation and dust accretion. Accretion as shown in Fig. 
[3] can only occur if the small debris created in catastrophic 
collisions is swept up by larger bodies rather than forming 
new planetesimals. Increasing tp by a factor of 10 still re- 
sults in accretion up to 80 km, but for a factor 100 growth 
stalls at 10 km. Similarly, the result depicted in Fig. |3] is 
robust to changes in td up to a factor of 10. 

We consider the a Centauri system in Figs. |4] and |5] 
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Figure 5. Evolution of maximum particle size (solid curve, top 
panel), mean particle size (dashed curve, top panel), planetesimal 
mass fraction (solid curve, bottom panel) and dust mass fraction 
(dashed curve, bottom panel) for the a Centauri system. Black 
curves are for a Centauri A, grey curves are for a Centauri B. 



accreting collisions. On the other hand, we have assumed 
dust accretion to be efficient. This efficiency depends on the 
size distribution that is produced in collisions and the ver- 
tical extent over which the fragments are distributed. Some 
of the fragments may be lost due to radial drift, which is 
ignored in the current model, but it may also help bringing 
more mass into the accretion friendly inner regions. 

We have worked in a 2D geometry for computational 
reasons. Although we have tried to scale the surface den- 
sity in such a way to get realistic collision time scales, it 
is important to realise that Tc,2d evolves in a different way 
with particle size than Tc. It is, therefore, not possible to 
have realistic collision time scales at all times. Furthermore, 
inclinations of the planetesimals may be excited by colli- 
sions, resulting in i = e/2, which decreases the collision time 
scale. However, since most collisions are destructive, it is 
likely that the remaining planetesimals will essentially orbit 
in the plane of the disc. Other important three-dimensional 



effects were pointed out by Xie & Zhou ( 2009 1 and Xie et al 



(20101, who showed that if the disc is inclined with respect 



to the binary plane, this may favour planetesimal accretion. 
Clearly, three-dimensional simulations are the way forward, 
and the results in this letter should be interpreted as a proof 
of principle, that planetesimal accretion might be possible 
in tight binary systems under the right circumstances. 



Results are shown for both binary components A and B, 
but since only the binary mass ratio changes between the 
two, the results are very similar. Again, we see growth up 
to 70 km in the inner region of the disc. Due to the fact 
that the system is more strongly perturbed than 7 Cephei, 
the accretion-friendly region has shifted inward compared 
to Fig. [2] The rough periodicity seen at early times in the 
bottom panel of Fig. [5] is due to different generations of 
planetesimals. The largest objects are formed inside 1.4 AU 
with significant growth inside habitable zone < 1 AU. 



6 DISCUSSION 

We present the first study of planet formation in binaries 
taking into account the physical effects of collisions. We 
show that these effects tend to favour growth towards large 
planetesimals in the perturbed system. Two main mecha- 
nisms have been identified. First, frequent collisions tend 
to prevent planetesimals from reaching their equilibrium 
orbits. If collision rates are high enough, this means that 
the accretion-hostile environment due to differential orbital 
phasing is never reached. Second, fragments produced by 
high- velocity collisions make up a large reservoir of material 
that is very easily reaccreted by the remaining planetesimals 
as they sweep through the disc on highly eccentric orbits. 

We have chosen parameters that in some respects 
should be unfavourable for planetesimal growth. Gas drag 
is strong throughout the disc in the models presented here. 
However, since gas drag plays only a minor role, changing the 
gas density to a more realistic power law in radius does not 
change the results. New planetesimals are formed on circular 
orbits and are, therefore, immediately capable of destroying 
larger bodies that are on eccentric orbits. Forming new plan- 
etesimals on eccentric equilibrium orbits would favour more 
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